Dispersion force for materials relevant for micro and nanodevices fabrication 
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The dispersion (van der Waals and Casimir) force between two semi-spaces are calculated using 
the Lifshitz theory for different materials relevant for micro and nanodevices fabrication, namely, 
gold, silicon, gallium arsenide, diamond and two types of diamond-like carbon (DLC), silicon carbide, 
silicon nitride and silicon dioxide. The calculations were performed using recent experimental optical 
data available in the literature, usually ranging from the far infrared up to the extreme ultraviolet 
bands of the electromagnetic spectrum. The results are presented in the form of a correction factor to 
the Casimir force predicted between perfect conductors, for the separation between the semi-spaces 
varying from 1 nm up to 1 /im. The relative importance of the contributions to the dispersion 
force of the optical properties in different spectral ranges is analyzed. The role of the temperature 
for semiconductors and insulators is also addressed. The results are meant to be useful for the 
estimation of the impact of the Casimir and van der Waals forces on the operational parameters of 
micro and nanodevices. 



I. INTRODUCTION 

In recent years, as the size of microelectromechanical 
systems (MEMS) decreased, the distance between mov- 
ing parts decreased correspondingly entering in the range 
of hundreds or even tens of nanometres. At such sepa- 
rations, the relatively large moving parts of MEMS ex- 
perience an attractive force due to quantum fluctuations 
of the electromagnetic vacuum that can be relevant for 
the determination of their operational parameters [l], Q • 
Borrowing a terminology frequently used in the realm of 
quantum field theory this attractive force is commonly 
referred to as Casimir force L 3j. Entering in the realm 
of nanoelectromechanical systems (NEMS), shorter sep- 
arations between moving parts are found which lie in 
the range of a few nanometres. In such case, it has be- 
come usual to consider the effects of the well known van 
der Waals force However, since the work of Lifshitz 
5] the van der Waals and Casimir forces can be treated 
in a unified form as resulting from the electromagnetic 
quantum vacuum fluctuations existing between two sep- 
arated dielectric bodies. The van der Waals (Casimir) 
force arises as a limiting case for short (large) separations 
between the bodies. Because in our analysis we are going 
to cover separations ranging from 1 nm up to 1 /im, we 
cover both limits. For that reason, and to avoid this some 
times confusing terminology, which includes the concept 
of retarded van der Waals force, we refer to the force here 
analyzed simply as dispersion force, a unifying terminol- 
ogy also found in the literature. 

The role that dispersion forces may play in the deter- 
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mination of the operational parameters of MEMS and 
NEMS has been analyzed theoretically in several in- 
stances [H, 0, E| and determined experimentally in, for 
example, references [3, 0, H[ . Different methods and ap- 
proximations were used in those analysis to determine 
the resulting dispersion force: pairwise summation of the 
interatomic Lennard- Jones potential Q , a rough distance 
independent correction factor to Casimir force between 
perfectly conduction plates [|| or the Lifshitz theory us- 
ing actual optical data P, Q for the determination of the 
force between silicon plates, or simple estimates based 
on the results for highly idealized cases In fact, be- 
cause the dispersion forces are highly dependent on the 
geometry of the interacting bodies and on their frequency 
dependent dielectric properties [IJ [l(| there is no sim- 
ple method of calculation that applies to every situation 
given reliable results. 

Nevertheless, for simple geometries, involving large 
surfaces with small deviations from the plane parallel 
geometry, a powerful, reliable and relatively simple ap- 
proach is to take into account the dielectric properties 
of the interacting bodies through the Lifshitz theory for 
the pressure resulting between two semi-spaces .and the 
geometry through Derjaguin approximation 0, 0, [l(| as 
explained in the next section. However, results based on 
the Lifshitz theory Q are available in the literature solely 
for three materials of relevance for micro and nanofabri- 
cation: gold and aluminium, for distances varying from 
less than 1 nm up to 1 /jm [l lj . and single crystal sili- 
con, for distances varying from 1 nm up to 1 /im Q. In 
the present article we expand this list including results 
for amorphous silicon, crystalline gallium arsenide, crys- 
talline diamond, two types of diamond-like carbon, cu- 
bic silicon carbide, amorphous silicon nitride, and amor- 
phous silicon dioxide, in the range 1 nm - 1 /an. 
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The article is organized as follows. In section [XT] we 
briefly summarize how Lifshitz theory is used in order 
to calculate the pressure between two dielectric semi- 
spaces, and how this result can be used to approximate 
the force between dielectric non-planar bodies. In sec- 
tion IIIII the results for the different materials we analyze 
are presented. We conclude in section [TV] 



II. 



LIFSHITZ THEORY AND DEJARGUIN 
APPROXIMATION 



According to the Lifshitz theory the dispersion force 
between two semi-spaces made of the same material and 
separated by a vacuum, is a function of the frequency 
dependent complex dielectric function of the material, 
e(u>). The resulting pressure at a separation d and at 
temperature T can be written as Q 
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were ks denotes the Bol tzmann consta nt, c the speed of 
light, K = K{i£) = VP 2 - 1 + e(iO, 6 = 27rk B Tl/h, 
and the prime near summation sign means that the ze- 
roth term is taken with the coefficient 1/2. The sum 
in equation |T]) [42| converges rapidly for temperatures 
above a few tens of Kelvin, however, for lower temper- 
atures the number of terms to be summed become too 
large and the use of the zero temperature approxima- 
tion becomes computationally more appealing. When 
fcfiT <C hc/d we can replace the sum by an integral and 
the pressure becomes, 
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At zero temperature and in the limit of perfectly con- 
ducting plates (e — * oo) equation ([2|) recovers the now 
classical result for the Casimir force acting between par- 
allel plates 
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while for dielectrics separated by short distances the pres- 
sure results to be 
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therefore recovering the van der Waals force (Ah denotes 
the Hamaker constant). 

Equations fl} and {2} depend on actual optical data 
through the complex dielectric function e(w) = £i(w) + 
ie2(oj). More specifically we calculate the dielectric per- 
mittivity along the imaginary axes appearing in both 
equations, e(i£), from e{oj) with the help of a Kramers- 
Kronig relation 
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In order to accurately calculate the pressure over a 
wide range of the separation d, it is necessary the knowl- 
edge of £2(0;) in a correspondingly wide range of frequen- 
cies, in such a way that the result of the approximate nu- 
merical integration performed to obtain e(i£) be as close 
as possible from the exact result. In order to fulfill this 
requirement we considered the optical properties of the 
different materials in a range of frequencies ranging from 
infrared (IR) up to extreme ultraviolet (EUV) or x-ray. 
For all materials we consider, there is data available cov- 
ering the interband transition region, in the visible and 
ultraviolet (UV) regions of the electromagnetic spectrum. 
This is the spectral region that gives the more important 
contribution to the dispersion force in the dielectric ma- 
terials we consider. However, for large separations the 
contributions from IR region have to be considered while 
at the shortest distances we consider contributions from 
EUV have to be taken into account, as is going to be 
evidenced by our analysis. 

For some materials we analyze there is tabulated opti- 
cal data covering the entire relevant spectral region which 
can be interpolated an numerically integrated. However, 
for silicon nitride and silicon dioxide we have to perform 
an approximation resorting to the Ninham-Parsegian rep- 
resentation [9| of e(i£) in the IR region. In this represen- 
tation e(i£) is calculated based solely on the absorption 
strength (Cj) and relaxation frequency (u>i) of the most 
relevant absorption peaks (in any spectral region, in prin- 
ciple) , 
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with the second term approximating the integral in equa- 
tion ||5j). Therefore, in calculating e(z£) for silicon nitride 
and silicon dioxide we add the contribution of each ab- 
sorption peak in the IR region. 

For gallium arsenide, silicon carbide, silicon nitride, 
and DLC, the optical data is limited to relatively low 
energies situated in the UV region of the spectrum. In- 
tending to improve the results for the shortest distances 
we consider we extend the data set by inclusion of the 
data on the index of refraction n and extinction coeffi- 
cient k provided by the Center of X-Ray Optics at the 
Lawrence Berkeley National Laboratory (T^. This data 
results from the compilation of both experimental and 
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theoretical data for elements with Z = 1 — 92 on the 
atomic scattering factors [l3| . The resulting n and k (we 
remind the reader that €2 = 2nk) are expected to be re- 
liable for energies above a few tens of electronvolts since 
at such energies the effects of atomic bounds in the solid 
state system, not taken into account by the atomic scat- 
tering factors, become sufficiently small. We checked this 
expectation by comparison with experimental data for n 
and k for several solids. The disagreement varied from 
less than 5% to crystalline silicon to more than 30% for 
Germanium, however, on the average the agreement was 
within 20% for energies above a few tens of electronvolts. 
We comment on the consequences of using such x-ray 
data in section Hvl 

The results for the pressure between semi-infinite par- 
allel planes can be used for the approximate calculation 
of the pressure between surfaces of different geometries. 
Whenever the relative inclination or the curvature of the 
surfaces are sufficiently small and the geometries are suf- 
ficiently simple the Derjaguin approximation can be used 
@> This approximation permits that the results 

presented in the next section be used for the calcula- 
tion of the dispersion forces for a large class of MEMS 
and NEMS made from moving parts other than parallel 
plates. 



III. RESULTS 

In this section we present the results for the disper- 
sion force. However, instead of the resulting pressure 
we present the finite conductivity correction factor, 77(d) 
which is a measure of the influence of the finite conduc- 
tivity on the Casimir force. This factor is defined as 
the ratio between the pressure given by the Lifshitz the- 
ory and the Casimir pressure between perfect conductors, 
equation It evidences more clearly the effects of the 
different optical properties for the different materials we 
analyze. 

For most materials we present 77(d) for different tem- 
peratures: K, 300 K and, 700 K. The chosen temper- 
atures are representative of cryogenic temperatures, am- 
bient temperature and the high temperatures that may 
be found inside integrated circuits, respectively. When 
only one temperature is considered we take T = K. 
Our results span the separation range 1 nm - 1/im. Wc 
do not consider larger distances because the dispersion 
forces become too small to be relevant in any practical 
MEMS or NEMS. On the other hand, the smallest sep- 
arations we consider are already found in sophisticated 
NEMS like the nanoresonator presented in [3], were a 
gap of only 20 nm exist between the moving parts. 

Besides the curve for 77(d), for the most relevant cases 
we also provide a least-squares polynomial fit to the data 
points that can be readily used. The fitting polynomial is 
of the form T](x) = a + Y2n=i a n xn , were a new variable 
is used instead of d, namely, x = log 10 (l / d) . The fitting 
for all cases are accurate within ~ 1% for d > 3nm, but 



present larger discrepancies for shorter distances, in some 
cases exceeding 4% at 1 nm. The order of the polynomial 
is the smallest one required to keep the accuracy ~ 1% 
for d > 3 nm, and varied from 6 to 7. 



A. Metals (Gold) 

Metals are widely used as a structural material for the 
fabrication of microsensors and microactuators, specially 
in RF MEMS, and were used in NEMS, for instance, on 
top of nanoresonators as conductive layers to provide a 
means of both induction and detection of motion [l5[ 
or to produce nanocantilevers [l6j |. among many other 
applications. Aluminium (Al) and gold (Au) are among 
the most widely used metals in micro and nanodeviecs, 
while nickel (Ni), cromium (Cr), titanium (Ti) and some 
metallic alloys were also used. 

Due to their high electrical conductivity, metals are 
the materials that best approximate the perfect conduc- 
tor boundary condition leading to the highest observable 
forces. For this reason they were used in all the recent 
experiments intended to measure the Casimir force with 
the highest accuracy [3[. The Casimir force for Al and 
Au was already investigated theoretically using the Lif- 
shitz theory [llj. We present the results for gold (Au) as 
representative of the metals used for micro and nanofab- 
rication to allow an immediate comparison to the results 
for the semiconductors and insulators we consider in the 
present work. Here wc use the same data set and further 
procedures as those employed in reference [llj obtaining 
the same results. The Casimir force was calculated at 
T = K using equation We do not consider higher 
temperatures because of the present controversy on how 
to specify the dielectric function for metals at zero fre- 
quency [l7| • However, for a good conductor it is expected 
that the correction due to the non-zero temperature ap- 
proaches that for perfect conductors, specially at large 
distances H, being only slightly greater. The tempera- 
ture corrected pressure between perfect conductors reads, 
to leading order in T, 
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with the separation d given in micrometres. According 
to this expression, at 300 K the force increases by just 
0.15% at a separation of 1 fim. At 700 K, however, the 
correction is considerable at the same distance, the force 
increasing by 4.6%. But it decreases very rapidly and at 
d = 500 nm the correction is of only 0.3%. 

In figure [T] we present £2 and 77 for Au. We note that 
77 ~ 1 at 1 /^m and it decreases rapidly for shorter separa- 
tions. Around 1 nm, 77 changes linearly in the log-linear 
graph indicating that the force varies with d 3 . This result 
is expected at short distances and characterizes the van 
der Waals force between two semi-spaces as seen from 
equation (fj|. The fitting function for gold at K is: 
r]{x) = 1.28706 xl0 3 -5.81088xl0 4 a; + l. 08899xl0 6 x 2 - 
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1.08370 x 10 7 x 3 + 6.03528 x 10 7 x 4 - 1.78202 x 10 8 a; 5 + 
2.17818 x 10 8 x 6 . 



B. Silicon 

Silicon in both its crystalline (c-Si) and polycrystalline 
forms is the most used material in MEMS and NEMS 
fabrication. Crystalline silicon is widely employed is pro- 
cesses involving bulk micromachining, while polysilicon 
in those where surface micromachining is required. Crys- 
talline silicon is the best understood semiconductor, and 
its optical properties were thoroughly investigated. High 
quality optical data is available from the far infrared up 
to x-ray energies. Here we use the tabulated data of ref- 
erence [3 for e 2 for energies ranging from 0.85 eV up 
to 3750 eV. In the IR region the interaction of ligth with 
the phonons is rather weak because silicon is a homopo- 
lar crystal and there are no active optical phonons. For 
this reason e 2 below 0.85 eV can be neglected [l8j . 

As for any polycrystallyne material, the precise opti- 
cal properties of polysilicon are highly dependent on the 
sample properties; different deposition methods and post 
deposition treatments result in different optical proper- 
ties [Til, and the optical constants for different samples 
are available only for the visible spectral region. How- 
ever, the available data on n and k indicates that the 
optical properties of any polysilicon sample are going to 
be between those of c-Si and that for amorphous silicon 
(a-Si) [l9[. Therefore the curves for e 2 and r\ for c-Si and 
a- Si are the two limiting cases for polysilicon. For that 
reason and due to the fact that optical data for a-Si is 
available in a wide energy range, we are going to present 
r\ for a-Si instead of specific polysilicon samples. We use 
the tabulated data for a-Si of reference [HI ranging from 
0.7 eV up to 48 eV. The contributions from the IR re- 
gion can be safely neglected as for c-Si and we do not 
extend e 2 for a-Si into the EUV region because our main 
interest is simply to have a reasonable bound for r\. As 
discussed further in sectionHVlthe lack of the optical data 
beyond a few tens of electronvolts do not significantly af- 
fect the results for silicon for separations larger than a 
few nanometres. 

In figure [2] we present e 2 and rj(d) for c-Si and a-Si. 
At T — K, rj for a-Si is greater than that for c-Si for 
all d, and we expected that ij for polysilicon samples be 
somewhere between the two curves. The effect of the 
temperature is rather small for c-Si up d ~ 300 nm and 
becomes significant for larger separations only at T = 700 
K. For a-Si, a similar dependence on T was observed. The 
fitting function for c-Si at K is r](x) = -1.72448 x 10 3 + 
7.63984 x 10 4 s- 1.40018 x 10 6 x 2 + 1.35873 x 10 7 x 3 - 
7.36383 x 10 7 x 4 + 2.11393 x 10 s a; 5 - 2.51214 x 10 8 x 6 
and for a-Si rj(x) = -1.76475 x 10 3 + 7.81535 x 10 4 a; - 
1.43181 x 10 6 x 2 + 1.38893 x 10 7 x 3 - 7.52482 x 10 7 x 4 + 
2.15945 x 10 8 x 5 - 2.56554 x 10 8 x 6 . 



C. Gallium Arsenide 

Gallium Arsenide (GaAs) has been widely used on the 
production of MEMS and NEMS despite its, in many 
aspects, unattractive mechanical properties due to its 
compensating attractive electrical and optical proper- 
ties, availability as comercial high quality single crystal 
wafers, and the relatively ease control of deposition and 
bulk micromachining, among other reasons [2(J. Only 
the crystalline GaAs (c-GaAs) has been used in MEMS 
and NEMS fabrication and its amorphous form is rarely 
cited as a potentially useful material. For these reasons 
and due to the lack of optical data for amorphous GaAs 
[l8T ] we analyze in this section solely the dispersion force 
for c-GaAs. 

Because c-GaAs has two different atoms per unit cell 
and has ionic bonds, it has active optical phonons in the 
IR region. Therefore, contrary to silicon, we cannot sim- 
ply disregard the optical properties in the IR region and 
we resort to the tabulated optical data collected in refer- 
ence [l8[ which spans the energy range from 0.01 eV up 
to 15 eV, which contains one IR resonance peak. Above 
15 eV, the data presented in [18| is only for the extinc- 
tion coefficient and do not include the index of refraction. 
However, experimental data on n and k above 9.5 eV is 
also available in 12 II] . Comparison between this data with 
that collected in [18j evidences considerable disagreement 
in the overlapping energy range corresponding to values 
of e 2 differing by more than 30%. The differences extend 
over larger energies when we compare solely k from both 
references. Due to this incompatibility we resort to the 
x-ray data to obtain approximate results for e 2 in the 
entire energy range above 15 eV. The x-ray data is avail- 
able above 30 eV and e 2 for the region between 15 eV 
and this energy was obtained by a third order interpola- 
tion including points below 15 eV and above 30 eV. We 
expect the error due to the use of interpolated data to be 
acceptably small since a relatively smooth variation of e 2 
in this region can be expected based on the data from 

In figure [3] we present e 2 and rj(d) for c-GaAs. We 
note how rj for c-Si and c-GaAs are quite similar with 
differences that rarely exceed 5%. The impact of the IR 
activity of GaAs on the dispersion force is rather small 
even for the largest distances we consider, leading to a 
dependence of rj on the temperature that is similar to 
that for silicon. For c-GaAs at K the fitting function 
is r)(x) = -1.45323 x 10 3 + 6.45764 x 10 4 a; - 1.18680 x 
10 6 cc 2 + 1.15457 x 10 7 x 3 - 6.27138 x 10 7 a; 4 + 1.80393 x 
10 8 x 5 - 2.14758 x 10 s x 6 . 



D. Diamond and DLC 

Diamond has been used on the fabrication of MEMS 
and NEMS in both its polycrystalline and amorphous 
forms (DLC). It is suited for devices operating in harsh 
enviroments due to its chemical inertness and high hard- 
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FIG. 2: Upper panel: £2 as a function of photon energy. Lower panel: r\ as a function of the separation distance d. Results are 
for silicon. 
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FIG. 3: Upper panel: £2 as a function of photon energy. Lower panel: i/asa function of the separation distance d. Results are 
for crystalline gallium arsenide. 



ness (2Cf and for the fabrication of high frequency micro 
and nanomechanical resonators due to its large Young 
modulus, of the order of 1000 GPa, and relatively low 
density (p ~ 3.5 g/cm 3 ). 

Because the optical properties of polycrystalline dia- 
mond are highly sample dependent and the optical data 
is restricted to small spectral ranges, as for polysilicon, 
we do not present an analysis for this form of diamond. 
Instead we consider crystalline diamond (c-diamond), 
whose optical properties were measured from far infrared 
up to x-ray energies, and that can serve as a reference ma- 
terial. We use the tabulated data for crystalline diamond 
of reference [l8[ ranging from 5.5 eV up to 800 eV. e 2 be- 
low 5.5 eV can be safely neglected because diamond is 
highly transparent in the optical region and as for silicon 
the contributions from the IR region are small duo to the 
absence of active optical phonons. 

DLC is a metastable form of amorphous carbon con- 
taining a significant fraction of sp 3 bonds 22]. This ter- 
minology encompasses several forms of amorphous car- 
bon including purely amorphous carbon (a-C) and its hy- 
drogenated alloys (a-C:H) and the sp 3 bonding rich tetra- 
hedral amorphous carbon (ta-C) and its hydrogenated 
alloys (ta-C:H) [22j ■ a-C and a-C:H have been widely 
used as protective coating for the disk and read/write 
head in magnetic hard disks and are currently being in- 
vestigated as potential materials for MEMS fabrication 
having its low temperature deposition as an important 



advantage over materials like SiC, polycrystalline and ul- 
trananocrystalline diamond and ta-C [22Ll23l]. Compared 
to other forms of DLC, ta-C is the one that best preserves 
the best qualities of diamond, possessing high mechani- 
cal hardness and chemical inertness. Due to this qualities 
ta-C is an extremely interesting prospective material for 
MEMS due to its superior wear-resistant qualities, resis- 
tance to stiction (i.e. a combination of stickiness and 
friction) and potential as a biocompatible material that 
could be used inside the human body for medical pur- 
poses without generating an allergic reaction 24 1 . 

The optical properties of the DLC can vary widely de- 
pending on the fraction of sp 3 bonds and hydrogen con- 
tent. We consider two representative DLC samples, one 
for ta-C with a large sp 3 content of 80% and another 
for a-C:H with a relatively low hydrogen concentration 
of 25%. The optical data for e 2 for the ta-C sample were 
extracted from reference [25| for the energy range 1.7-40 
eV. For energies above 40 eV we used x-ray data calcu- 
lated considering a sample made from pure carbon and 
having the same density as measured in [25[ p = 3.0 
g/cm 3 . Because there is a mismatch between the experi- 
mental and x-ray data in the range 30-40 eV of the order 
of 15% we actually use the x-ray data for energies above 
50 eV softening the mismatch over the 10 eV gap between 
40 and 50 eV using numerical interpolation. The opti- 
cal data for a-C:H were extracted from reference [26| and 
spanned the energy range between 0.7 eV up to 30 eV. 
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As for ta-C we observed a considerable mismatch at 30 
eV between experimental and x-ray data (the x-ray data 
was calculated for a sample having four carbon atoms for 
one hydrogen and the density obtained experimentally of 
p = 1.75 g/cm 3 .) and used a 10 eV gap in order to soften 
the mismatch. Therefore, x-ray data was used above 40 
eV. With regard to the IR activity of ta-C and a-C:H, it 
was experimentally verified [2jJ that both forms of DLC 
have low absportion in the IR region, and therefore £2 in 
the IR can be neglected. 

In figure [5] we present £ 2 and 77(d) for c-diamond and 
the two forms of DLC. The resulting correction factors 
r] present an interesting behavior. At short distances, 
below approximately 100 nm the Casimir force for c- 
diamond is slightly larger than that for ta-C, however 
for separations above 100 nm an inversion occurs. This 
inversion can be easily understood. The force at short 
distances is influenced mostly by the strength of the res- 
onances in the UV region which are clearly weaker for 
ta-C as compared to c-diamond. At larger distances the 
force is influenced by lower frequency modes for which €2 
is larger in the case of ta-C . It is interesting to note also 
the decrease in the dispersion force for a-C:H relative to 
c-diamond and ta-C, making a-C:H an interesting choice 
to minimize the Casimir force. 

The fitting function for c-diamond at K is the 7th 
order polynomial T)(x) = -6.88550 x 10 3 + 3.53738 x 
10 5 x - 7.74634 x 10 6 x 2 + 9.37393 x 10 7 a; 3 - 6.77091 x 
10 8 x 4 + 2.91985 x 10 9 x 5 - 6.96200 x 10 9 x 6 + 7.08214 x 
10 9 x 7 . 



phonon. We use the data collected in [l8| between 0.01 
eV up to 0.4 eV which covers the IR region containing 
the resonance just mentioned. For energies above 0.4 eV 
and up to approximatelly 20 eV we use the data on £2 of 
reference [29[ . Above 30 eV we resort to x-ray data, and 
use an interpolating function in the 20-30 eV range. 

In figure [5] we present £2 and r)(d) for 3C-SiC. In or- 
der to determine the impact of the IR activity of the SiC 
on 77(d) we present in figure \E\ rj calculated at K with- 
out the contributions from the IR region (indicated by 
the label "No IR" in the figure). The result is consider- 
ably smaller than that for the full spectrum at K for 
distances above approximately 100 nm, demonstrating 
the importance of the IR activity for this material. More 
specifically, the relative increase due to the IR activity on 
the force is about 21% at a separation of 1 71m. For non- 
zero temperature we observe a larger relative increase of 
X] as compared to the materials analyzed previously, spe- 
cially for distances near 1 /im. were it increases by 25% 
at 700 K compared to K. This larger increase can be 
attributed to the contribution of the IR region. The dif- 
ferent consequences on 77 of the IR activity for GaAs and 
SiC can be attributed to the much stronger IR resonance 
peak of the later (by about a factor of six) and to its 
occurrence at a higher energy around 0.1 eV, compared 
to 0.033 eV for GaAs. 

As for c-diamond the fitting function for 3C-SiC at 
K is a 7th order polynomial, namely, -q{x) = —4.67473 x 
10 3 + 2. 36272 x 10 5 a; - 5.08848 x 10 6 x 2 +6. 05401 x 10 7 a; 3 - 
4.29817 x 10 8 x 4 + 1.82144 x 10 9 x 5 - 4.26710 x 10 9 x 6 + 
4.26431 x 10 9 x 7 . 



E. Silicon Carbide 

Silicon carbide (SiC) is a material with very attrac- 
tive properties for MEMS and NEMS applications sim- 
ilar to those of diamond. Its mechanical strength, high 
thermal conductivity, ability to operate at high temper- 
atures and extreme chemical inertness, make SiC attrac- 
tive for MEMS and NEMS applications in harsh enviro- 
ments, both as structural material and as coating layer, 
and for high frequency micro and nanomechanical res- 
onators. SiC has been used on the fabrication of MEMS 
and NEMS in both its crystalline and polycrystalline 
(poly-SiC) forms [2(3] and the potential usage of its amor- 
phous form has been investigated (28|. In its crystalline 
form, SiC is a polymorphic material that exist in cubic, 
hexagonal, and rhombohedral polytypes. However, only 
for the cubic polytype (3C-SiC) there exist adequately 
measured optical data [181 ] in a sufficiently wide range 
of energies ranging from 0.01 eV up to 20 eV [HI, [29| . 
Therefore, due to the lack of data for polycrystalline and 
amorphous SiC, as well as for the non-cubic polytypes 
we restrict our analysis and consider the Casimir force 
for crystalline 3C-SiC. 

Because SiC is a heteropolar semiconductor it is active 
in the IR region. 3C-SiC presents a strong IR resonance 
around 0.0987 eV due to an active transversal optical 



F. Silicon Nitride 

Silicon nitride (SiN) is widely used in MEMS for elec- 
trical isolation, surface passivation, etch masking, and as 
a mechanical material [20], being in the last case also 
used in NEMS, for instance, on the fabrication of nano- 
electromechanical oscillators 30]. Silicon nitride is usu- 
ally employed in its amorphous form and frequently it is 
also non-stoichiometric and may contain significant con- 
centrations of hydrogen. Good optical data is available 
in a wide energy range for the amorphous stoichiometric 
silicon nitride (a-SiaN^ [3lJ, and for crystalline stoichio- 
metric (3— SiaN4 j32|. Here we focus on the results for a- 
SiaN4 because SiN is usually employed in its amorphous 
form in MEMS and NEMS. 

For a-SisN4 we use the data from [3l[ which spans the 
energy range 4-24 eV. There are two resonances in the IR 
region whose effects we introduce using equation ([6]) with 
Ci = 1.08 and w x = 0.9 x 10 14 rad s"\ and C 2 = 2.37 
and cl>2 = 1-64 x 10 14 rad s _1 [HI]. For energies above 30 
eV we use x-ray data calculated for stoichiometric SiaN4 
with density p — 3.2 g cm -3 , and use cubic interpolation 
between 24 and 30 eV. We note that the £2 generated 
using x-ray data has very good agreement with the ex- 
perimental data for /3-Si3N4. 
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FIG. 4: Upper panel: £2 as a function of photon energy. Lower panel: r/ as a function of the separation distance d. Results are 
for diamond. 
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FIG. 5: Upper panel: £2 as a function of photon energy. Lower panel: r/ as a function of the separation distance d. Results are 
for silicon carbide. 
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In figure [6] we present £2 and 77(d) for a-SiaN4. From 
the results for 77 it can be concluded that the resulting 
Casimir force is smaller than that for the previously an- 
alyzed materials. The influence of the IR activity is also 
noticeable, as for SiC, having a larger relative importance 
on 77, that corresponds to a 34% increase on the force at 
1 /zm and K. The effect of temperature is also large for 
a-SiaN4, leading to a 33% increase on the force at 1 /im 
and 700 K as compared to K. 

The fitting function for a-SiaN4 at K is a 7th order 
polynomial, namely, 77(2) = -4.57954 x 10 3 + 2.33430 x 
10 5 x - 5.07148 x 10 6 x 2 + 6.08838 x 10 7 a; 3 - 4.36263 x 

10 8 x 4 + 1.86622 x 10 9 x 5 - 4.41385 x 10 9 x 6 + 4.45362 x 

10 9 x 7 . 



G. Silicon Dioxide 

Silicon dioxide, in its amorphous form (a-Si02), is usu- 
ally known for its applications in microelectronics as an 
electric insulator. More recently it has been widely used 
in the fabrication of MEMS and NEMS as a sacrifi- 
cial layer material and has also found application as a 
structural material, for instance, for the fabrication of 
large suspended membranes and grating light valves [34| . 
Quartz, the crystalline form of SiC>2, has also been used 
as a structural material for the fabrication of MEMS, 
specially in microfluidics. 

For both a-SiC>2 and quartz there is good optical data 
available in a wide ran ge o f energies and we resort to the 
most recent works. In [35j the data covers the range 1.5- 
42 eV for both forms of SiC>2, while in f3|| the energy 
range 65-3000 eV is covered for the amorphous form. 
Due to the availability of recent data in a wider range 
for a-SiC>2 and to its almost omnipresence in MEMS and 
NEMS we analyze the Casimir force for this form of sil- 
icon dioxide. However, due to the similarity between 
the optical properties of crystalline and amorphous Si02 
evidenced in [35j, our results are expected to be approx- 
imately valid for quartz. 

Complementing the experimental data from [35l[36l| we 
use x-ray data between 42 eV and 65 eV, which has shown 
a very good agreement with both data sets between 30 
eV and approximately 100 eV where data superposition 
exists. From the IR region we considered the existence of 
three resonances whose effects we introduce using equa- 
tion 



with d = 0.829 and uji = 0.867 x 10 14 rad 



C* 2 = 0.095 and w 2 = 1-508 x 10 14 rad s" 



and 



C* 3 = 0.798 and u> 3 = 2.026 x 10 14 rad s" 1 [33j. 

In figure [7] we present £2 and rj(d) for a-Si02- From 
the results for 77 it can be concluded that the resulting 
dispersion force is considerably smaller than that for the 
previously analyzed materials, a result attributable to 
the largest band gap of about 8.3 eV [37]. The influence 
of the IR activity is also the strongest amongst all the 
materials we analyzed increasing the force by 74% at 1 
fxm and K. The effect of temperature is also rather 
large, leading to a 52% increase on the force at 1 /im and 



700 K as compared to K. 

The fitting function for SiC>2 at K is a 7th order 
polynomial, namely, rj(x) = -1.84016 x 10 3 + 9.41606 x 
10 4 x - 2.05374 x 10 6 x 2 + 2.47524 x 10 7 x 3 - 1.78058 x 
10 8 x 4 + 7.64622 x 10 8 x 5 - 1.81524 x 10 9 x 6 + 1.83827 x 
10 9 x 7 . 



H. Mixed materials 

As our last result, we present 77 for combinations of 
different materials in each semi-space. In this case we 
use a generalization of equation (|2|) for the pressure be- 
tween semi-spaces made from different materials found, 
for instance, in (38|. As representative of all possi- 
ble combinations we consider the following three com- 
binations of Au, c-Si and a-Si02: Au-Si, Au-SiC>2 and 
Si-SiC>2. They represent the three possible combina- 
tions conductor-semiconductor, conductor-insulator, and 
semiconductor-insulator. The results for 77 are presented 
in the figure [5J where, as a reference, we also reproduce 
the results for Au, Si and Si02- It can be seen that, as 
could be expected, the resulting force is always between 
that produced for each material separately. 



IV. DISCUSSIONS AND CONCLUSION 

In performing our calculations we used the most re- 
cent and complete optical data available for each mate- 
rial in order to derive the most accurate results. How- 
ever, even the best tabulated optical data relies on results 
from different experiments based on different samples and 
methodologies covering different, and in many cases, not 
overlapping spectral re gion s that must be interpolated or 
modelled theoretically jl8l . l3l| . Besides there may exist 
a significant disagreement between experiments, of the 
order of tens of percent, as pointed in section IIII CI for 
GaAs and as verified by the authors also for the case of 
3C-SiC and a-Si02 when the results from references [l8[ 
and [1^] and [35[ and [3!| were compared, respectively. 
Therefore an analysis of the expected uncertainties in our 
theoretical predictions is necessary. 

The result for gold must be considered taking into ac- 
count the thorough analysis performed in [4(| that has 
shown that great uncertainties, around 5%, must be as- 
sumed for the predicted force between gold plates. How- 
ever, the analysis performed for gold can not be extended 
for semiconductors and insulators. We have firstly to 
rely on the analysis performed by Aspnes and Studna 
(4l| in 1983 in order to state that for the well studied 
c-Si, GaAs, c-diamond and a-Si02 the uncertainty on 
the latest optical data for £2 is about 2%. This is an 
estimate based on the analysis of the convergence of dif- 
ferent data sets, and means that better measurements 
should not result in changes for the optical data larger 
than this uncertainty. This is the only uncertainty esti- 
mate we can perform since, as a rule, the optical data 
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FIG. 6: Upper panel: £2 as a function of photon energy. Lower panel: r\ as a function of the separation distance d. Results are 
for amorphous silicon nitride. 




FIG. 7: Upper panel: £2 as a function of photon energy. Lower panel: r\ as a function of the separation distance d. Results are 
for amorphous silicon dioxide. 
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TABLE I: Maximum uncertainty on r\ due to variations on 
£2- Approximate values for small and large separation d. 
Material(<5e2) Sr/ (Small d) 8-q (Large d) 



c-S (2%) 1.6% 


0.9% 


c-GaAs (2%) 1.7% 


0.9% 


c-diamond (2%) 2.0% 


1.4% 


3C-S1C (5%) 5.3% 


3.2% 


a-Si 3 N 4 (5%) 5.8% 


3.6% 


a-SiO-2 (2%) 2.9% 


2.0% 



presented in the references we used are not accompained 
by an analysis of theoretical or experimental uncertain- 
ties. For a-Si, 3C-SiC and a-Si3N 4 larger uncertainties 
on £2 should be assumed because few experiments were 
performed to measure their optical properties. We sug- 
gest that an uncertainty of at least 5% be assumed. In 
the case of DLC, the dependence of optical properties on 
sample preparation make it worth to consider our results 
only as a good estimate to the dispersion force. 

The resulting uncertainty on r\ (Srj) due to the un- 
certainty on £2 (^2) depends on the separation d and 
material under consideration. For instance, for silicon Srj 
is always smaller than <5e 2 , being within 40% to 80% of 
its value. However, for 3C-SiC and a-SiaN4 Sr] results to 
be larger than be-i for short distances, reaching a maxi- 
mum of 8r] — 5.8% for <5f2 = 5% in the case of a-SisN4. 
For larger distances Srj decreases to approximately 70% 
of (5e2 for both materials. As a reference we present in 
table U the maximum and minimum Sij as a function of 
Se2 for the most relevant semiconductors and insulators 
we considered. The general trend for the uncertainty is 
that it decreases rapidly from its maximum at 1 nm to a 
value close to its minimum for d ~ 100 nm and is almost 
a constant for larger d reaching a minimum at 1 /im. 

For c-GaAs, 3C-SiC, a-SiaN4 and DLC another source 
of uncertainty results from the use of x-ray data. As al- 
ready mentioned in sectionUTlcomparison between the ex- 
perimental and x-ray data for several materials, like alu- 
minium, copper, silicon and germanium, evidences that 
x-ray data is precise within about 20% for energies above 
30 eV. The actual disagreement could be smaller, since 
the x-ray data tends to agree well (within 10%) with the 
most recent, and presumably more precise, experimental 
data for materials like silicon and silicon dioxide. We 
can assume conservatively that Se2 for such energies is 
about 20%. In order to estimate how this uncertainty 
affects i], we firstly note that for the semiconductors and 
insulators we are considering 62 above 30 eV amounts to 
at most 6% of i] for d ~ 1 nm, 3% at d ~ 10 nm and 
contributes less than 1.5% at 100 nm. Therefore, the 
extra uncertainty due to the x-ray data is at most ap- 
proximately 1.2% at 1 nm, 0.6% at 10 nm and 0.3% at 
100 nm, and can be neglected for larger separations. It 
has to be noted, that this extra uncertainty results in an 
error that is much smaller than the error resulting from 
the neglect of €2 above 30 eV, therefore justifying the use 
of the x-ray data when no experimental data is available. 



The above uncertainty analysis demonstrates that a 
complex relationship exists between the uncertainty in 
the optical data and the theoretical prediction of the 
Casimir force. Therefore, for precise experiments a care- 
ful analysis must be performed for each material and ex- 
perimental condition. However, the results we present 
are useful for theoretical estimates of the impact of dis- 
persion force on the operational parameters of MEMS 
and NEMS. 

We also addressed the role of temperature on the ex- 
pected dispersion force since MEMS and NEMS operat- 
ing in real conditions may be subject to high tempera- 
tures that alter the force. For all materials but SiC-2 the 
effect of the ambient temperature of T — 300 K on the 
force is rather small. At T = 700 K for silicon, gallium 
arsenide, and diamond the effect of temperature is not 
so severe but exceeds the approximate result for metals 
by at least a factor of two for d ~ 1/xm. However, for the 
materials with strong IR activity SiC, S13N4 and Si02 
the effects of the temperature are already noticiable at 
shorter distances for T = 700 K and are rather large for 
d in the range of hundreds of nanometres. At this point, 
it is worth to mention that no experiment performed to 
date has measured the dependence of the dispersion force 
on temperature but, according to our results, materials 
with strong IR activity are the most adequate for such 
an experiment due to the greater sensitivity of the force 
on temperature. 

Our analysis also evidences how significant can be the 
contribution for the dispersion force coming from the IR 
region through the results for SiC, S13N4 and Si02- The 
effects of the IR activity are noticiable already at dis- 
tances in the range of a few tens of nanometres and in- 
crease the dispersion force by tens of percent at larger 
separations even at K. Therefore, due to its potential 
effects on the dispersion force, knowledge of the IR opti- 
cal behavior of a given material is mandatory when this 
force is to be calculated for distances ranging from tens 
up to hundreds of nanometres. The IR spectral region 
can be neglected only after a careful analysis. 

When the dispersion forces are to be calculated for 
short distances, smaller than about 100 nm, knowledge 
of the optical properties in the UV region is required. 
We considered tipically the optical properties up to a 
few hundreds of electronvolts, corresponding to the EUV 
spectral region. However, as expected for semiconductors 
and insulators, the more important contribution to the 
dispersion force over the entire interval considered for the 
separation d comes from the interband transition region, 
usually extending from ~ 1 eV up to ~ 20 eV. 

To conclude, we note that the dispersion forces for the 
different materials relevant for micro and nanodevices 
fabrication can vary considerably. When large disper- 
sion forces are desirable, good conductors or small band 
gap semiconductors should be used. However, in order 
to minimize the dispersion force, insulators are recom- 
mended. The results for diamond and DLC are an ex- 
ample on how the dispersion force can be taylored by 
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d(m) 



FIG. 8: rj as a function of the separation distance d. Results are for two different materials on each semi-space. The continuous 
lines are (downward) for Au, c-Si and a-SiC>2. 



chemical or physical changes of a given material as a Brazil), 
strategy to control the dispersion forces in MEMS and 
NEMS. 
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